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ABSTRACT 

Wc propose a new model for the description of ultra-short flares from TeV blazars by compact 
magnetized condensations (blobs), produced when red giant stars cross the jet close to the central 
black hole. Our study includes a simple dynamical model for the evolution of the envelope lost by 
the star in the jet, and its high energy nonthermal emission through different leptonic and hadronic 
radiation mechanisms. We show that the fragmented envelope of the star can be accelerated to Lorentz 
factors up to 100 and radiate effectively the available energy in gamma-rays predominantly through 
proton synchrotron radiation or external inverse Compton scattering of electrons. The model can 
readily explain the minute-scale TeV flares on top of longer (typical time-scales of days) gamma-ray 
variability as observed from the blazar PKS 2155—304. In the framework of the proposed scenario, the 
key parameters of the source are robustly constrained. In the case of proton synchrotron origin of the 
emission a mass of the central black hole of Mbh ~ IO^Mq, a total jet power of Lj « 2 x 10^^ erg s""'^ 
and a Doppler factor, of the gamma-ray emitting blobs, of ^ > 40 arc required. Whilst for the external 



inverse Compton model, parameters of Mbh ~ 10 M©, Lj w 10 
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Subject headings: galaxies: jets — 
(PKS 2155-304) 



Gamma rays: galaxies 



erg s ^ and the 6 > 150 are required. 
- BL Lacertae objects: individual 



1. INTRODUCTION 

The flux variability of very high energy (VHE) gamma- 
rays on minute timescales detected from the BL Lac ob- 
ject PKS 2155-304 (lAharonian et al.| |2007| ) and Mkr 



501 (Albert et al. 20071 ) challenges the standard scenarios 
suggested for the e xplanation of the nonthermal proper 



ties o f TeV blazars ( [Begehnan et al. 2008; Giannios et al 
^009 ) . The extremely short duration of the hares impose 



severe constraints on the size of the gamma-ray produc- 
ing region, of 



/' < cr' ~ 3 X lO^V^ cm. 



(1) 



where I' and Tj = r'/lO'^ s are the proper production 
size and the variability timescale in the frame of the jet 
respectively and c is the light speed. The proper vari- 
ability time-scale, r', is connected to the variability in 
the observer frame, r, by the relation 



r 



(2) 



where 6 is the Doppler factor of the moving source (the 
blob): 

^ " r(l-/3cos(a)) ■ 

Here the bulk Lorentz factor, P, accounts for the rela- 
tivistic transformation of time, and (1 — /3cos(a)) is re- 
sponsible for the kinematic shrinking of the duration of 
the radiation and f3 = v/c. 

Relativistic jets ejected from the central engines are 
common phenomena for different types of active galac- 
tic nuclei (AGN). In particular, apparent superluminal 



speeds /3app = /3sina/(l — (3 cos a) (in units of the speed 
of light c) as high as ~ 40 have been detected for radio 
com ponents on (projected) scales of ^ 1 — 10 pc f see 



e.g. [lorstad et aL| |200l| , |2005| ; [Marschcr et al.| ^OlOj ) 



blazars - AGN with jets directed to the observer. This 
implies very large Lorenz factors of bulk motion given 
that P > /3app. 

Because of the large bulk Lorenz factors of TeV blazars, 
the condition Equation (0) allows significant relaxation 
of the r equirement supplied by Equa tion (|l|). In par- 
ticular, Lcvinson &: Bromberg| (2008) argued that if a 
perturbation is produced by the central engine, its size 
should exceed the gravitational radius of the black hole 
in the observer frame. Consequently, the proper size 
of the production region is expected to be larger than 
Prg, where rg = GMbh/c^ ~ 1.5 x 10^^AfBH,8 cm is 
the gravitational radius of a black hole (BH) of mass 
A/bh.8 = Mbh/^O^Mq. In this case, the variability time- 



scale T2 = r/lO s imposes a strict upper limit on the 
gravitational radius: 



r„ < 3 X 10 T9— cm 

8 "P 



(4) 



Thus, the detection of gamma-ray flux variability r ~ 
200 s constrains the BH mass, to be Mbh, 8 < 1. In 
reality, since the main energy release occurs in the inner 
parts of the accretion disk of radius r ~ lO^g, the mass 
of the black hole should be close to IO^'Mq. Generally, 
in the case of an extreme Kerr black hole, the energy 
release takes place at the gravitational radius. However, 
even in this case one needs an entire rotation period for 
an effective energy release, i.e. the characteristic time 
cannot be much shorter than 27rrg/c. This implies that 
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even in this case the upper bound of the BH mass of 
IO'^Mq cannot be significantly relaxed. 

This conclusion is true, in particular, for the model of 
internal shocks. Note, however, that it is based on the as- 
sumption that the perturbation (the reason of the flare) 
originates in the central engine. Therefore it cannot be 
unco nditionallv extended to other possible scenarios as 



is cla imed by Dermer et al. ( 2009). Indeed, Equation ( 4) 
is not valid if perturbations are produced by an exter- 
nal source, e.g. by plasma condensations (often called 
"blobs") which do not have a direct link to the central 
black hole. Such blobs can be produced, in particular, 
by intera ctions of stars with the base of the jet, as pro- 
posed by Barkov ct al. (2010) to e xplain the TeV flare s 



2010), 



of M 87 on scales of days (see also Araudo et al 
where the interaction of gas clouds from broad line re- 
gions are discussed). The jet power of M 87 is relatively 
10*4 



modest, L 



The results of Barkov et al 



.J _ -Lv. ergs ^ ^_ 

show that, while this power is sufficient to blow- 

up the envelope (atmosphere) of the star which initially 
has been pulled out by the tidal force of the BH, such 
a jet appears to be not sufficiently powerful for acceler- 
ation of the gas cloud to relativistic velocities. Actually 
this works in a positive direction for M 87, given the 
large aspect angle of the jet. Otherwise, the gamma-ray 
flux could not be observed because of the Doppler de- 
boosting effect. On the other hand, the suggested mech- 
anism of formation of hadronic blobs in the jets cannot 
apply to powerful gamma-ray blazars unless the blobs 
are accelerated to Lorentz factors F > 10. Remarkably, 
this can be realized in a quite natural way in powerful 
jets with I/j > 10^^ erg s""'^. Interestingly, such power- 
ful jets can ablate the star atmosphere without help of 
tidal forces (the interesting implications of this effect are 
discussed below). Moreover, the powerful jets can drag 
and disrupt the star's envelope into an ensemble of blobs 
moving with large Lorentz factors which, from the point 
of view of explanation of very short time variability is a 
quite comfortable situation. 

Another important aspect of the short time variabil- 
ity is related to the efficiency of acceleration and ra- 
diation mechanisms. Currently the most conventional 
approach for modeling of VHE emission production in 
active galaxies is based on the inverse Compton (IC) 
scattering of relativistic electrons, the soft target pho- 
ton field being either the synchrotron radiation of same 
electrons (the so-called Synchrotron-self Compton (SSC) 
model), or provided by external sources (EIC model). 
The apparent advantage of IC models is the combina- 
tion of two factors: (1) the acceleration of electrons to 
relatively modest energies (< 1 TeV) can be effectively 
realized within different acceleration scenarios; and (2) 
these electrons radiate readily in interactions of ambient 
radiation and magnetic fields. Nevertheless, while the IC 
models allow rather satisfactory explanations of the en- 
ergy spectra and variability patterns of many blazars in 
general, the parameters used to fit some specific objects 
appear incompatible with the parameters defined from 
observations. Moreover, the observed short variability 
time-scale demands conditions which appear to be quite 
uncomfortable, in terms of the strength of magnetic field 
and related consequences concerning the strong deviation 
from equipartition between the energy density of rela- 



tivistic electrons and magnetic fields. The requirement 
of weak, (generally less than 1 G) magnetic fields is one 
of the key postulates of the IC paradigm of gamma-ray 
production in blazars. Moreover, in the case of some ob- 
jects with unusually hard source spectra (after correction 
for intergalactic absorption), such as lES 0229-f 2 00, the 



magnetic field is required to be as small as ImG ( Tavec 



chio et al. 2009). The magnetic field in the blazar jets 



can be reduced to such small values only at very large 
distances from the central engine, namely > 10^* cm. 
Although this idea has some observational support re- 
l ated to the transpar ency of blobs in the radio band 
(Marscher et al, 2008), it is likely that regions of highly 
variable gamma-ray emission ar e located at smal l er dis - 
tances from the central engine dTavecchio et al.||201C| ). 
In particular, the EIC models require the location of any 
gamma-ray emitter to be located closer to the BH into 
the so-called Broad Line Regions (BLR), i.e typically at 
distances R ^ 10^^ — 10^^ cm. This implies that any 
IC model can be realized only if one finds a way dra- 
matically reduce the magnetic field in the jet. Although 
this cannot be excluded (e.g. because of reconnection 
of th e B-field (Komissarov et al. ^007a ; Giannios et al 



2009 



tive 



3ulk acceleration of pl asiiiaT Komissarov et al 



McKinney fc Uzdensk>-| pOiq) or due to the effee^ 



Tchckhovskoy ct al.|poTo| )), strong magnetic fiek 



201C; 



s ex- 



ceeding 1 G remain a more favored option, as long as 
we deal with strong jets on sub-parsec scales. In this 
regard, the models which invoke high energy protons for 
the production of gamma-rays passes certain advantages 
despite a quite p opular view that they are not effective 



emitters (sec e.g. Sikora 2010). Actually this is true only 



for proton-proton and proton-photon interactions. What 
concerns the synchrotron radiation of protons, with a key 
assumption on the acceleration of particles with a rate 
close to t^^ ^ ceB/E, where E is the proton energy, 
coupled with a strong magnetic field between 10 to 100 
Gauss, and a large Doppler factor, (S > 10, it is that can 
provide relevant acceleration and radiation timescales, 
as well as explain the extension of gam ma-ray spectra to 
TeV energies (sec e.g. Aharonian 2000| ). 



2. BLOBS IN RELATIVISTIC JETS 

2.1. AGN Jet - Red Giant interaction 

Below we discuss the distinct features of the interac- 
tion of red giants with AGN jets in the specific case of 
powerful blazars, see sketch of the scenario in Figure |l|. 
Originally, the scenario of AGN Jet - red gian t f RG ) in- 



Barkov et al 



teraction has been suggested by ^ 

the explanation of VHE observations of M87 - a non- 
blazar type nearby AGN with a large jet viewing angle 
of > 20° and a modest jet power, of Lj^t — 10'*'* erg s^^. 
It was demonstrated that if disturbed by tidal forces RG 
penetrates the jet, the ram pressure of the jet in M87 
would be sufficient to remove the outer layer of the RG. 
This leads to the formation of a dense cloud within the 
jet which, in combination with effective particle accelera- 
tion, can trigger gamma-ray production through proton- 
proton interactions which passes secure of the required 
spectral and temporal properties. This model allows de- 
tectable gamma-ray fluxes because of the proximity of 
the source and the non-relativistic speed of the blobs 
(otherwise the radiation from M87 would be de-boosted 
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given the large aspeet angle of the jet). 

The AGN jet - red giant interaction (JRGI) has very 
specific and important features in the case of powerful 
blazars, where the ram pressure can be as high as 



10'^ L j _4qz^^ 9 dyncm 



(5) 



Here ij,46 = ij/lO^^ergs^^, zn = z/lO^^cm and 9-i = 
0/0.1 are the jet power, the distance from the BH, and 
the jet opening angle, respectively. Such a high ram 
pressure can blow-up the outer layers of the stellar at- 
mosphere, even from a non-disturbed RG. The mass of 
the removed layer can be roughly estimated from the 
balance of the jet ram force to the gravitational force, 
Pre.ra'^Rl ^ GAMM^/Rl. Here M, and are the RG 
mass and radius, respectively. This gives the mass of the 
cloud stripped from the star by the jet: 



AM 



6 X IQ^^L^aqz^^OzIR 



*,2M~U^ (6) 



where R^^2 = R^/W^Rq and i\/*,o = MJMq are the RG 
parameters expressed through the solar mass (Mq) and 
radius (i?©). This estimate illustrates the feasibility of 
formation of a cloud due to the JRGI process. The next 
important issue is the acceleration of the cloud in the 
jet. In what follows, we assume that initially the ejected 
matter from the RG envelope forms a (quasi) spherical 
cloud of radius comparable to the RG size, Tc — R* ■ 

To be accelerated, the matter residing in the jet first 
should be heated, thus the cloud expands enhancing the 
interaction between the cloud and the jet. A significant 
expansion occurs on a time-scale of f pvp = Atr-r.. where A 



is a constant of order o f 5 (see e.g. Gregori et al 



al 



Pittard et al.| |20iq ), a"iid 



2000 



ire. 



Nakamura et alj E006 

2rc/cs is the sound crossing time. Here Cg is the sound 
speed in the shocked cloud, which can be estimated as 
Cs « 7g-Pram/Pc ~ 7gPram47rr3/(3Afc) (7g = 4/3 is the 
plasma adiabatic coefficient). This gives the following 
estimate for the expansion time: 



t, 



cxp 



At,, = 5 X lO^i?, j/^Af^'/'gij ]fzi70-i s, (7) 

where M,^26 = M,/I0^^g is the normalized mass of the 
cloud. The orbital velocity of a star around the black 
hole, 



orb 



GA/bh 



4 X 10^ A/, 



ri/2 -1/2 

BH,8^17 



cm s 



results in a crossing time of the jet of 



BH 



= 3 X lQ^z%'^i 



n r-1/2 
8 S- 



(8) 



(9) 



The condition tcxp < ^jc defines an upper limit on the 
mass of the cloud, which can expand in the jet: 

Afc,sc 



<3 X 1029ij,46i?*,22l7A^Bff,8g- 



(10) 



If this constraint is fulfilled, then at the stage of hydro- 
dynamical expansion, the cloud can incre ase its linear 



size b y a few orders of magnitude (see e.g. [Pittard et al 
^OlOD whilst remaining in the jet 



The cloud will be trapped by the jet if it would be 
sufficiently accelerated along the jet axis. Thus, the con- 
finement condition has the following form: Vz = a^tjc > 



K)rb/£ where az — PramT^r^/Mc. Using Equations (^) 
and (m one obtains the cloud capture condition: 



Afcjc < 3 X 10''Li,i6rl,,M^i, g 



(11) 



where rc.15 ~ rc/lO^^cm is the size of the cloud after the 
hvdrodynamical expansion. Importantly, Equations (|lO| ) 
and (|ll]) provide upper limits, which exceed significantly 
the expected mass of the blown-up layer given by Equa- 
tion (0). Thus, even if the ablated stellar matter remains 
as one blob, this cloud will still be trapped in the jet and 
accelerated up to sub-relativistic velocities (i.e. ~ 0.1c). 
In fact, independent of the initial conditions, during the 
acceleration phase the cloud is expected to be c rushed 
into hundreds of small blobs ( Pittard et al. 2010| ). This 
relaxes significantly the above obtained conditions, since 
the ensemble of blobs can be more easily picked up by 
the jet flow. 

2.2. Relativistic Stage 

At the relativistic stage, the dynamics of the cloud is 
described by the following equation: 



(12) 



where T, is the Lorentz factor of the cloud, lj z6 
is cylindrical radius of the jet (see Appendix ^ for the 
derivation of Equation (|l2|)). Let us introduce the follow- 
ing notations g = Pc/Fj and y = z/zq (and dy « cdt/zo). 
Here zq is the distance from the BH to the point where 
the RG penetrates the jet, i.e., we adopt the initial con- 
dition 5 <C 1 at z = Zq. Several simplifications, in partic- 
ular, assuming that Pj = const on the cloud acceleration 
time scale, and dt = dyzo / c, allow presentation of Equa- 
tion ([l2[) in the form: 




where 



dg 
dy 

D 



1 



D 



— ^9 — ' 

9 J y 



46l2p3zoc3A^c 



(13) 



(14) 



Equation (|T^) allows an analytical solution, which defines 
y as a function of g and D. 

In this paper we do not specify the origin of relativistic 
particles in the cloud. However, assuming that particle 
acceleration is a result of a strong interaction between 
the cloud and the jet, we may conclude that high Lorentz 
factors do not support the production of non-thermal ra- 



diation. Indeed, as it follows from Equation ([L3D, the 
cloud-jet interaction intensity decreases with cloud ac- 
celeration. The apparent intensity of the non-thermal 
phenomena can be roughly described by the luminos- 
ity correction function Fe-, which accounts both for the 
Doppler boosting, i.e. F, cx P^, and for the interac- 
tion intensity, i.e. F, cx — /z^- For the sake 

of clarity, we consider the correction function in dimen- 
sionless form F, == q^(l/q^ — g'^)/y'^. Since the solu- 
tion of Equation (|l3) relates g and v, F, is, in fact, a 
function of one variable. In Figure we show it as a 
function of g (left panel) and as a function of the vari- 
able >f = t(l - V,)/to = 2D{y - l)p2(l - V,) (right 
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panel), corresponding roughly to the observation time 

(here Vc = vl — TV^ is the blob velocity) . The values of 
the parameter D = 0.1, 1, 10 and 100 are used in both 
panels of Figure p. It can be seen that in the case of 
D > 1, the solution depends weakly on the parameter 
D, and the non-thermal activity of the blob is expected 
to have a narrow peak of duration: 



to 



c D 2r2 ■ 



(15) 



The maximum of the correction function occurs at 
ffmax ~ 0.8. 

In the case of D <^ 1, the situation is quite different. 
Namely, the expected nonthermal activity has no pro- 
nounced peak, thus such a blob cannot produce a high- 
amplitude flare. In this regard, we can formulate the 
condition D > 1 as a requirement for a flaring episode. 
This condition can be reformulated as a upper limit on 
the cloud mass: 



(16) 



The above relation depends not only on the properties 
of the cloud (its size and mass) but also on the jet power 
and Lorentz factor. Thus, for quantitative calculations, 
one needs detailed information about the dynamics and 
properties of the blazar jet. 

Although the process of jet formation is not fully un- 
derstood, recent hydrodynamical studies of differ ent sci- 
entif ic groups s how that the Blandford-Znaicck ( Ruffini 
& W tlson||l975t |Lovelacc|[l976| ; iBlandford fc Znajel|Tfrt ) 
process may be at work in AGN, and suggest a concept 
of magnetically accelerated jets. Thus, the jet base is 
expected to be strongly magn etized and likely magnet- 
ically dominated at z < 1 p c ( Komissarov ct al. ^007b 



Barkov fc Komissarov| |2008| ; peskin| p010| ). This imme 
diately gives the magnetic held strength of the jet in 
laboratory frame 



Si 



1/2 



1204/46^1"/' 



'Ii'G- (17) 



During the jet propagation, the magnetic field energy can 
be transformed to the bulk kinetic energy. At the linear 
s tage a simple relat ion defines the bulk Lorentz factor 
(iBeskin fc Nokhrina] p006|) 



UJ 

4r„ 



(18) 



Finally, th e opening angle of the j et is expected to be 

Combining Equa- 



e ^ 1/r 

tions dl 



( Kornissarov ct al, 
and 



meie 
^009| ). 



, one can estimate the magnetic field 
in the jet comoving frame. 



Br 



2 fU 



1/2 



12Z;^/l/j G . 



(19) 



In Figure ^, the typical magnetic field and bulk Lorentz 
factors of the jet are shown for three different distances 
from the BH. These values are in good agreement with 
observed values of magnetic field on parsec scales in 



AGNs ( 


Lobanovl 


199^; 


Savolainen et al. 


200s ; 0' Sullivan 


fc Gi 


ibuzda 


2009 


; Sokolovsky ct al. 


201( 





Using Equation (^8|) , one can present Equation ( p^ ) in 
the form 

Afc,.c < 0.5 X 102'^Lj,46r2 igFrf ^A/Bj^ ^ g. (20) 

Note that at this late stage, the cloud can be already 
significantly expanded with a radius of rc,i5 3> 1. The 
extreme value of rc can be achieved at Tc ~ w: 



Merc < 2 X 102^Lj,46MBH,8rj:i.5g- 



(21) 



This upper limit is more robust than the constraints 
given by Equations (|6l),(0), and ([11]), we note however, 
that it concerns the mass of the blob, but not the mass 
of the ablated stellar atmosphere. 

2.3. Energy Budget of the Cloud 

Below we consider the general requirements to the sce- 
nario in the context of blob's radiation efficiency. These 
constraints have a quite basic character and are not re- 
lated to a specific radiation mechanism. Obviously, the 
gamma-ray production mechanisms impose additional 
requirements, concerning e.g. the density of relevant tar- 
gets in the form of gas, radiation or magnetic field. We 
discuss the impact of specific radiation mechanisms in 
Section ^ Here we try to find a generic link between 
properties of a blob, which is responsible for the non- 
thermal emission, to the parameters of the AGN, i.e. the 
mass of the central engine and the jet power. 

Emission detected from blazars arc significantly en- 
hanced due to Doppler boosting. The apparent {L^) 
and intrinsic (Lsc) luminosities are connected through 
the well known relation: = LscS^ (here is Doppler 
factor of the production region). On the other hand, the 
intrinsic luminosity of the blob can be expressed as a 
fraction ^ of the power transferred by the jet to the blob: 



(22) 



The parameter ^ ^ 1 accounts for the overall efficiency 
of of transformation of the absorbed jet energy to non- 
thermal emission through acceleration and radiation of 
relativistic particles (note that the righthand-sidc of the 
Equation (|22| ) refers to the quantities in the observer 
frame, whilst the lefthand-side corresponds to the blob's 
reference frame). For small aspect angles, S « 2rc, we 
thus obtain the following simple relation: 




(23) 



which has a few interesting implications. In particular, 
one can estimate the size of the blob: 



1/2 



(24) 



Given the standard bulk Lorentz factors, Fj sa a;/(4rg), 
and the maximum of correction function, max (Fe) « 0.4, 
together with conventional normalizations, one obtains 



r inl47ir rl/2 J--1/2, — 1/2 

5 X 10'*MBH,8-^^,47^j,46 C-1 



(25) 



Here we use quite a high normalization for the efficiency, 
of^_i=^/0.1 
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Another important estimate can be obtained for the 
maximum apparent luminosity of the blob. It is achieved 
when the blob eclipses the whole jet, i.e. rc ~ to. In this 
case, one obtains: 



2x 10*Viij,46rf,i.5ergs-i. 

-,1.5 



(26) 



where Fj^i s — Fj/lO . Note that the apparent non- 
thermal luminosity of the blob is proportional to T?. 
Remarkably, even for relatively modest values of the jet 
Lorentz factor, Fj ~ 10, and the conversion efficiency, 
^ ^ 0.01, i-y is comparable to the jet power. 

2.4. Time variability 

The fast variability of TeV gamma-ray emission of 
blazars, which can be as short as a few minutes as re- 



ported for PKS 2155-304JAharonian et al.| |2007| ) and 
Mkr 501 ( Albert et al. 2007 ), is a key observational fact 
which should be addressed by any model of TeV blazars. 
In addition to a general (standard) statement about the 
strongly Doppler boosted gamma-ray emission produced 
in very compact regions close to the central black hole, 
any dedicated model should provide intrinsic reasons for 
variability (characterizing the scenario as a whole) and 
offer radiation mechanisms with adequate cooling times. 
If the cooling time-scale appears to be too long, an alter- 
native source of the variability is required. For example, 
a change of the production site velocity may lead to a 
strong change in the apparent luminosity. Indeed, since 
the production region is to be Doppler boosted, a varia- 
tion of the Doppler factor may be a plausible reason for 
strong variation of the observed flux. 

In this section, we discuss the variability scales related 
to this effect, and its implication to the JRGI scenario. 
We note that the relevant time-scales do not depend on 
the cooling time of the emitting particles. In particular, 
the obtained size of the blob in Equation (GS) and the 
jet Lorentz factor provide following lower limit on the 
variability time-scale: 



^ > _^ « 4 X lO'z-^''L^^%L7U'CrM^'l, s , (27) 



which appears to be close to the observed one, and can 
be significantly shorter in the case of powerful jets. 

2.4.1. Duration of the blob-jet interaction 

The principal variability scale in the JRGI scenario is 
related to the duration of the effective interaction of the 
cloud with t he j et and is determined by the function 
(see Section 2.2). Since the model requires very effective 
acceleration of particles, with a > 10% efficiency of the 
energy transformation to nonthermal particles, the shape 
of the function, Fo, (see Figure ^ can be treated as the 
time profile of particle acceleration with a characteristic 
timescale: 

2^rV^ . (28) 



At 



Note that in the extreme case, when the blob eclipses the 
entire jet, i.e. Zgfl^/r^ ~ 1, the characteristic time-scale. 
At, depends only on the jet Lorentz factor Fj and power 
I/j, as well as on the mass of the cloud Mc'. 



The total apparent energy of electromagnetic radiation 
which can be emitted by the cloud can be estimated from 
Equations (|2l) and (poh: 



1050e-iMc,25ni.5erg. 



(30) 



For the values expected in this scenario (which were 
also used for normalization of Fj, Lj and Mc in Equa- 
tion (p9|)), the cloud-jet interaction can be quite brief; 
- shorter than the detected variability of the ultrafast 
flares of PKS 2155-304 and Mkn 501. Moreover, for 
small mass clouds, it can be as short as 1 sec. Obviously, 
this time-scale corresponds to the flare rising interval, 
while fast emission decay requires short radiative or adi- 
abatic cooling of the emitting particles or rapid changes 
in the blob's Doppler factor. 

2.4.2. Helical structure of relativistic jet 

Generally, in powerful jets, the change of the Doppler 
factor is unavoidable. Indeed, since the matter in the 
jet moves along the dominant magnetic field lines, which 
are expected to be helical, the velocity of the blob should 
have both poloidal and azimuthal components . The az 
imuthal velocity can be as high as Vd, cr i c/uj ( Deskiii & 



2007b| [Beskii^ |201CU ; 



Nokhrina 2006; Komissarov et al 

where ric ^ 4?'^ is light cylinder radius, therefore, during 
the motion, the angle towards the direction to observer 
can be changed. We discuss the corresponding variabil- 
ity pattern in Appendix where it is shown that this 
effect (which reminiscent to some extent, the operation 
of a 'revolver'), can lead to a change in the flux by a fac- 
tor of 2 if the blob has ~ 0.5c and turns around the 
jet axis by an angle of 7r/4. In the case of magnetically 
driven jet using Equation (|l^) we can get « c/Fj or 
if Fj > 3 this mechanism cannot to explain variability of 
blazars. 

2.4.3. Collision of blobs 

The interaction of the jet with a massive cloud leads 
unavoidably to the formation of a large number of small 
blobs, which may gain an additional chaotic velocity 
component. Let us assume that the chaotic velocity is 
comparable to the sound speed in relativistic gas, i.e. 
Vs ~ c/-\/3, which is larger than c/2, that is enough to 
explain strong radiation variability (see for details Ap- 
pendix Due to interactions with each other, the blobs 
can change their speeds on timescales tg ^ 2rc/vs leading 
to the variability on a timescale 



T « tJ5 > 5 X 10^Mbr^sL\%L7U^c1^^62^ : 



(31) 



At«60Fj,i.5Lj;46^/c.25S 



(29) 



where Doppler factor is determined as S2 = 

Equation (^) shows that the interaction time-scale 
in the JRGI scenario could be very short allowing, in 
the ease of comparably short particle cooling time, flar- 
ing episodes of duration ~ 100 s. If the particle cool- 
ing mechanism cannot provide the required energy loss 
rate, the nonthermal flux variability can be caused by a 
change in the production region Lorentz factor, see e.g. 
Equation (|3l]). However, we note that in this case the 
variability on timescales as short as 100 s would require 
a rather specific combination of several principal param- 
eters. Thus, a radiative mechanism with short cooling 
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time remains still a quite feasible requirement for mod- 
els intended to explain the fast variability observed in 
blazars. 

In the following section we discuss the efficiency and 
features of major radiation mechanisms related to both 
protons and electrons. Although the jet composition is 
still debated, the conventional approach attributes the 
non-thermal activity of AGN to a lepton IC mechanism. 
In the suggested scenario even if the primary content of 
the jet is leptonic, the ablated cloud itself may provide 
protons for the acceleration process. In this paper we do 
not discuss the specific mechanisms of particle accelera- 
tion, but simply assume that both electrons and protons 
are effectively accelerated during the interaction of the 
blob with the jet. 

3. RADIATION MECHANISMS 

In this section we discuss the applicability of differ- 
ent radiation mechanisms responsible for the gamma-ray 
emission of blazars in the context of the shortest variabil- 
ity timescale of order of 100 s observed during strongest 
flares of PKS 2155—304. The spectral energy distribu- 
tion (SED) of the source has a typical shape for such 
objects with two pronounced humps. The low energy 
component peaks in the optical-UV-soft X-ray band at 
V 10^^ Hz, whilst the high energy bump has a max- 
imum in the VHE band. The location of the gamma- 
ray maximum is measured only in the low state of the 
source, during simultaneous Fermi-HESS observations in 
2009, which revealed a broad maximum between 10 GeV 



and 100 GeV ( [Aharonian ct al.| |2009b| ). Although the 
exact position of the gamma-ray maximum is not yet 
measured in the high state of the source, the o bserved 



spectral flattenin g during the July 2006 flares ( Aharo 



nian et al. 2009a) indicates a tendency of the extension 
of the region of the flux maximum, but most likely not far 
beyond 100 GeV. Note that for the distance to the source 
of = 540Mpc {z = 0.116), the attenuation of gamma- 
rays of energy <C 1 TeV in the extragalactic background 
light (EBL) is negligible. Therefore throughout this pa- 
per we will assume that the gamma-ray peak in the SED 
of this source is located around 100 GeV. The average 
apparent gamma-ray luminosity of the 2006 July flares 
was at the level of ^ 10^^ ergs^^. During the giant July 
28 flare, the source was not monitored in the X-ray en- 
ergy band, but the simultaneous observations of the next 
night, also characterized by strong flares, conducted with 
the H.E.S.S., Chandra and Bronberg optical telescope, 
revealed that the luminosity of the object in optical, UV 
and X-ray energy bands was an order of magnitude lower 
compared to the gamma-ray luminosity. These rather 
general properties of the SED, which include peak loca- 
tion and flux ratio, allow us to derive some important 
constraints on the production mechanisms. 

3.1. Gamma rays associated with electrons 
3.1.1. SSC model 

In the SSC models the high energy gamma-rays are 
produced by relativistic electrons through IC scattering 
of synchrotron radiation of the same electron population. 
Generally, this model satisfactorily explains the basic fea- 
tures of gamma-ray blazars. However, the ultrafast flares 
of PKS 2155—304 pose severe constraints on the param- 



eters characterizing the gamma-ray production region. 
Generally, the IC scattering proceeds in the Thomson 
when hvE^ <C rn^c 5'^, where 5^1. Thus, in 



regime, wJien hvhj^ <C m c' 
the co-moving reference frame. 



(32) 



- 4 X 10-1V67', 

where 7 is the electron Lorentz factor, iJ-y.ii = 
£;^/100GeV and viq = i^/lO^'^Hz are the peak energies 
of IC and synchrotron components in the observer ref- 
erence frame. Then, the electron Lorentz factor can be 
estimated as 



7 



5 X lO'' 



J^16 



1/2 



(33) 



The strength of the magnetic field in the co-moving 
frame can be defined from the location of the synchrotron 
peak, 

= 6 X Hz, (34) 

where Bq = B/IG. Here, for order-of-magnitude esti- 
mates, we adopt that the maximum of the vF^, distri- 
bution of synchrotron photons occurs at energy 1.33wcf|, 
where 

ieBE^ 
2 m'^c° 

is the synchrotron characteristic frequency, while the 
maximum of the distribution is located at a lower 
energy of 0.29a;c. Thus, one obtains: 



Bn 



0.7 ly 



16^7,11^ 



(35) 



The ratio of the IC and synchrotron peak luminosities, 
/, is another important parameter characterizing two- 
hump SEDs: 

/ = i«L^^. (36) 



^syn 



Wb 



In the SSC model, taking into account the constraint 
on the size of the production region imposed by the ob- 
served variability of the time-scale in Equations (|l|-||), a 
lower limit for the co-moving energy density of the target 
photons can be obtained: 



Wph > 10 



erg cm 



(37) 



where Lx,46 = Lx/lO^^ erg s~^ is the apparent syn- 
chrotron luminosity. For the given luminosity ratio of 
the observed high and low energy components, one finds 



900 



1/4 



(38) 



Thus, for SEDs, typical for ultrafast flares of 
PKS 2155—304, the Doppler factor of the relativistically 
moving gamma-ray source should exceed S ~ 500. Note 
that this condition is stronger than the constraint on 
the Doppler boosting imposed by the condition of the 



janmia-ray transparency of the source (Begelman et al 
2008|). 



^ This energy is in fact used in Equation for IC 
in the Thomson regime 



scattering 
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The ordcr-of-magnitude estimates of Equations ( |35| ) 
and (IH); obtained for SSC scenario, suggest that short 
flares from PKS 2155—304 should be produced at large 
distances from the BH (given the weak magnetic field 
and large Doppler boosting factors). This requirement 
is, in fact, very constraining for JRGI scenario, since the 
jet ram pressure in this region appears to be extremely 
small, viz. 



Pr 



ram, SSC 



5x 10-3i/4g^-2^dyn,cm-2, (39) 



which is not enough (by far) to ablate the required 
amount of stellar material. 

The above estimates show that SSC models meet se- 
vere limitations in the framework of JRGI scenario due 
to the required weak magnetic field. We note that the 
above severe constraints are basically due to the very 
small magnetic field. In principle, one can assume that 
the magnetic field in the blob is much weaker then in the 
jet, which could improve the effect of the SSC mecha- 
nism. However, in any case. Equation ( |38|) can be sat- 
isfied only at very large distances, where the jet ram 
pressure is very small and to oblate and accelerate the 
stellar atmosphere is difficult. 

3.1.2. Model of external photon field 

The main difference between EIC and SSC models is 
that in the former one the gamma radiation is dominated 
by the scattering of electrons on low-energy photons of 
external origin. Obviously, in such a case the energy 
density of the external photon field should exceed signif- 
icantly the energy density of the synchrotron photons, 
i.e. the photon energy density in the jet vicinity should 
fulfill the following requirement 



Wcxt > 10 



10 



LxA6 _3 

;m erg cm , 



(40) 



which follows immediately from Equation (|37|). In fact, 
even a very weak external photon field can fulfill this 
requirement, given a very strong dependence on the jet 
Lorentz factor. Another limitation on the external pho- 
ton field can be derived from the ratio of the IC and 
synchrotron peaks (see Equation (|3^)) which together 
with Equations (|l^) and (|l9|) give the following luminos- 
ity limit: 



LcKt > 3 X 10*^/Lj^46Zi/MBH,8ergs" 



(41) 



which is likely available in the vicinity of powerful 
blazars. Since this photon field may remain unde- 
tectable, the energy of the target photon is, to a large 
extent, a free parameter. Thus, Equation (32) is not 
valid in the EIC case, and the model parameters are 
less constrained than in the SSC scenario. In particular, 
this allows us to relax the requirement of weak magnetic 
field, which is crucial for SSC models in the framework 
of JRGI scenario. To make quantitative estimates, one 
needs to assume some basic properties of the external 
field, namely the typical photon energy e r-y ~ e/1 eV 
and its luminosity L^xt- Then Equations ( ^2[ ) and (|35| ) 
obtain the following form: 



7,11 



10 eeV7 <J : 



(42) 



and 



(43) 



Together with Equations (|l^ and (p^, one can solve 
these equations and represent 7 and Cov a-s 



6eV « 104/46^7,11^4h.8^iV^17'^' > (44) 



and 



7 ~ 1.0 X lU ly^g Z-^-j JWbH,8-'^7,46 



(45) 



An important characteristic of the model is the emitting 
particle cooling time in the jet reference frame: 



Coi = 3x lOVW^'s, 



(46) 



where wp = wb + Wph = (/ + 1)wb is the comoving frame 
energy density of the target fields in units of erg cm~^ . 
Equations (|l^) and (^) result in the following cooling 
time: 



Col = 3 X 10^(1 + /) 



-17/4.-3/4 , .-1/4 -1/2^ 
'^17 -^j,46 ^"BH,8 "^16 



(47) 



This fast cooling allows us to relate the observed vari- 
ability on timescale of several hundred seconds, with the 
relativistic particle radiation cooling. We note as well 
that given the relatively low energies of the electrons, 
£^0 1 — 10 GeV, and the strong magnetic field [B ~ 
few G), the acceleration of these electrons can be easily 
realized. 

Finally, we have to note that the obtained values 
should fulfill the condition of IC scattering in the Thom- 
son regime: 76, 
e.g. Equation 

only). This yields the following requirement for the in- 
teraction point: 



2 (since some used relations, 
are valid in the Thomson regime 



Zi7 » h-^^QlVl, 



p4/3 -2/3 
BH, 8-^7, 11*^16 



(48) 



In the case that this condition is not true, the interac- 
tion of electron with target photons occurs in the Klein- 
Nishina regime and a few additional effects, such as 
Klein-Nishina electron losses and gamma-gamma atten- 
uation, have to be taken into account. Actually, since in 
this case, the Klein-Nishina losses have to be the domi- 
nant cooling mechanism, the gamma-gamma absorption 
is unavoidably large. Indeed, the Klein-Nishina cooling 
time is 



KN 



E 

Ekn 



1 



t^KN 



lOVsEj s. 



(49) 



where rig^t and ckn are the target photon density in the 
jet frame and Klein-Nishina cross-section, respectively. 
On the other hand, t he gamma-garn ma optical depth can 
be estimated as (see Derishev 2009| ) 



'77 



a^~f « 40AfBH,8T2 ^ , 



(50) 



where next and cr^^ are the target photon density in the 
laboratory frame and pair-production cross-section. To 
derive Equation (|50|), we have used the approximate re- 
lation (Tjj w 2crKN and Equations (|l8|) and (^7|). Thus, 
it is rather unlikely that short flares can be produced in 
the Klein-Nishina regime on the external photon fleld. 
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3.2. Gamma rays associated with protons 
3.2.1. pp and P7 interaction 

The production of gamma-rays from interactions of rel- 
ativistic protons with the surrounding gas is one of the 
major processes in high energy astrophysics. This pro- 
cess is effective in relatively dense environments, namely 
when the pp cooling time, tpp = s does not ex- 

ceed other characteristic times (here n is the target pro- 
ton density in units of cm~^). In the case of the blob in 
the jet, the most relevant times are the interaction time 
of the blob with the jet given by Equation (|2^ ) (which 
can be treated as the acceleration time of protons), and 
the escape time of protons from the blob. The number 
densitv of protons in the blob is constrained by Equa- 
tions(|l|) and (|5|). For the parameter values expected 
in the scenario, the number density of protons does not 
exceed 10^ cm~'^. The corresponding pp cooling time 
of > 10^^ s is too long, and does not leave any room 
for the explanation of the variability of TeV radiation 
on any observed time-scale. Even the assumption that 
the variability is caused by other reasons, e.g. due to the 
adiabatic cooling or change of the Doppler factor, cannot 
help much since in this case, an extremely low efficiency 
represents an unavoidable argument against this process. 

The efficiency of gamma-ray production can be much 
higher through other channels related to interactions 
with the radiation and magnetic fields. In both cases 
the gamma-ray production rate increases dramatically 
with the energy of protons, and achieves reasonably high 
efficiency if the protons are accelerated to energies of 
10^^ eV or beyond. For the compact blobs with lin- 
ear dimensions severely constrained by the hour-long or 
shorter variability timescales, the energy of protons can 
achieve such high energies only when the particle acceler- 
ation proceeds (1) at a rate close to the theoretical limit, 
iacc ^ Ti^/c and (2) in the presence of a magnetic field 
as large as 100 G (Aharonian ^000 ). Whilst for these 
conditions both the acceleration and synchrotron cool- 
ing times in the frame of the blob can be as short as Ih, 
the fast cooling of protons via photomeson interactions 
require very dense photon fields at mm and far-infrared 
wavelengths. On the other hand, the density of the radi- 
ation field (of internal or external origin) is constrained 
by the condition of transparency of the production re- 
gion for the VHE gamma-rays, implying that the optical 
depth regarding the photon-photon pair production can- 
not significantly exceed unity. This condition, coupled 
with the condition of a large magnetic field amplitude 
of _B ^ 100 G, makes the cooling time of protons via 
photomeson processes significantly lo nger compared t o 
the proton synchrotron cooling time (Aharonian 2000). 
Although, formally one can "construct" a model with ex- 
treme parameters, where the photomeson processes could 
compete with the proton synchrotron cooling, below we 
will focus our treatment on the production of gamma- 
rays via synchrotron radiation. 

3.2.2. Proton synchrotron radiation 

Protons of extremely high energy and large magnetic 
fields strengths are the two conditions which make proton 
synchrotron an effective radiation mechanism. If these 
conditions are satisfied, the spectrum of synchrotron ra- 
diation can extend to the gamma-ray domain with a char- 



acteristic energy (Aharonian 2000| ) 



E. 



'7,11 



IB2E 



19 ' 



(51) 



where Eig = E'/IO eV is the proton energy, and B2 = 
B/lOO G is the strength of the magnetic field. The po- 
sition of the peak depends strongly on the maximum 
energy of protons, which is determined by the balance 
between the particle acceleration and cooling rates. It is 
convenient to present the acceleration time of the pro- 
tons, independent of the specific mechanism of accelera- 
tion, in the form: 



t. 



i]{E)ri^ 



lO^EigB^'r^iE) s , 



(52) 



where tl — E/eB = 3 x 10-'^''i?2£'i9 cm is the so-called 
gyro-factor. We note that this value is remarkably close 
to the required size of the blob in the JRGI scenario (see 
Equation ( [25 1 )). The dimensionless parameter, ri{E) > 1, 
in Equation {\)2\) characterizes the acceleration efficiency; 
the most efficient acceleration occurs for t] = 1. It is be- 
lieved that the accelera tion by relativistic shocks (see e.g. 
Achterberg et al. pOOl ) or at the ann ihilation of the mag- 
Ha 



netic field lines ( Haswell et al. 1992 ) occurs in the regime 
when r] ^ 1. In a more general context, the relativis- 
tic outflows found in AGN and GRBs, which presently 
are considered as most feasible sites of the production of 
the highest energy cosmi c rays, can accelerate p rotons to 
10^° eV only if - 1 ( [Aharonian et al ||2002| ). In the 
comoving frame, the characteristic synchrotron cooling 
time of the protons is: 

ATTm'^C^ . „ -, 

W = ^r-z:^ = 3 X 10^ B^'E^,' s. (53) 



aTUi^EB^ 



From comparison of the acceleration and cooling rates 
one finds a maximum energy of protons of; 

-Bp,,nax = 1.7 X 10i'S2~'/'7?-i/' cV , (54) 

and a corresponding maximum energy of synchrotron 
photons of; 



i^7,max = 4OO77-I GeV 



(55) 



If synchrotron radiation is emitted by a blob with a 
Doppler-factor, 6, the peak of the synchrotron radiation 
is shifted to 



E^ 



400?7"M GeV 



(56) 



One can see that for a Doppler factor oi 6 < 100, syn- 
chrotron radiation can extend to TeV energies provided 
that the proton acceleration rate is close to 77 < 100. 

Another important constraint can be obtained through 
so-called Hillas criterion, i.e. a requirement for the size 
of the acceleration site to be larger than the gyroradius of 
the highest energy particles. Substituting Equations (pj| ) 
and (p4) into the Larmor radius, tl ~ E/eBc, one ob- 
tains 



3/2 j- 



-1/2 r -1/4^-1/2 fl/2 



17 -^7,47 ^3 



1 ^^BH,8 < 0-1 



(57) 



Here the size of the acceleration site was assumed to be 
defined by Equation (|25|). A significantly more severe 
constraint for the proton synchrotron models comes from 
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the cooling time requirement. Namely, the expected vari- 
ability time scale can be estimated from Equations ([l^), 
TqI), (|5|) and (H): 



''psyn ~ ~ 2 X 10 IJi 



4„l/2^.1/2 



BH 



46 



S . 



(58) 



Thus, in the case of powerful jets, ij > lO'^^ergs^^, the 
proton synchrotron mechanism can guarantee variability 
on time-scales of several hundred seconds. 

Although a detailed study of broad-band SED in the 
framework of proton synchrotron scenario remains out- 
side of the scope of this paper, below we outline some ex- 
pected features. In particular, one may expect that the 
acceleration of protons is accompanied by the accelera- 
tion of electrons, whose particle population may create 
a detectable nonthermal emission component. Due to 
more effective synchrotron losses, the maximum energy 
of electrons is rather small: 



= ^^^Vc^^^ TeV 



(59) 

The corresponding synchrotron peak appears at an en- 
ergy, which is by a factor of m-p/rrie smaller than the 
peak energy of proton synchrotron (see Equation (p6|)): 



IbQ-q^^ MeV. 



(60) 



Thus, in this scenario, the synchrotron peak of electrons, 
additionally boosted by the value of the Doppler factor, 
5 ~ 50, is expected at GeV energies, but not in the tradi- 
tional X-ray radiation band. In this regard the so-called 
synchrotron peak, which in TeV blazars is located in the 
UV to X-ray band, needs a new (non-standard) inter- 
pretation in the framework of our model. Finally, it is 
important to note that the flux ratio of proton to electron 
synchrotron peaks depends on the ratio of the injection 
rates of those particles, which is an highly uncertain pa- 
rameter. 

The peak of soft radiation in the region of 10^^ Hz can 
be produced by secondary electrons produced in inter- 
actions of VHE gamma-rays with soft photons. A mod- 
est absorption of VHE gamma-rays cannot be excluded. 
Moreover, in the case of blazars with extremely hard TeV 
source spectra (i.e. after correction for the intergalactic 
absorption), the energy- dependent absorption can be a 
natural explanation of the unus u ally hard gamma-ray 
spect ra (Aharonian et al. 2008; ^acharopoulou et al 

poTTI) 



The energy of the absorbed gamma-ray photon is 
shared between an pair, so each electron on aver- 
age gets the half of the original energy, i.e., E^ ^ E^/2. 
Generally, synchrotron radiation of secondary electrons 
produced inside the jet will be Doppler boosted, there- 
fore this amplified component of synchrotron radiation 
can be detected by the observer. The peak energy of the 
secondary synchrotron radiation is quite sensitive to the 
Doppler factor of the blob and the energy spectrum of 
the parent gamma-rays. For the standard values used 
in this paper, the peak of the synchrotron radiation can 
varv in a quite broad interval b etween UV and hard X- 
ray ( Zacharopoulou et al. [2011 ). 

Thus within the suggested scenario, we expect a broad- 
band SED consisting of three synchrotron peaks of differ- 
ent origin located at the keV, MeV/GeV and TeV bands. 
Schematically, such a SED is shown in Figure ^. While 



the amplitudes of the MeV and TeV synchrotron peaks 
are determined by the total energy accelerated in the 
form of electrons and protons, respectively, the intensity 
of the low energy synchrotron peak is determined by the 
fraction of absorbed TeV (proton-synchrotron) gamma- 
rays. 

4. APPLICATION TO PKS 2155-304 

In the case of the July 2006 flares of PKS 2155-304 
the total energy of the nonthermal radiation detected 
during the burst was about i?tot ~ L^I^t « 10^^ erg. 
According to Equation (pG), such an energy release can 
be produced by an ensemble of clouds satisfying to the 
following condition: 



j4-5 



10. 



(61) 



This requirement can be easily fulfilled given the large 
mass of material, that can be accelerated by the iet. In- 
deed, for the mass obtained through Equation (pi|), the 
above condition is reduced to: 



ij,46 > 0.5 



1 



BH,i 



(62) 



- j,1.5 



The detected short variability of r ^ 200 s implies 
some additional constraints on the system parameters. 
The causality condition constrains the size of the blob 
Tc ~ 5tc ~ Fj Tc (here we assume the size of the produc- 
tion region to be 2rc and 5 = 2Fj). This allows one to 



relate the variabilit:v 
through Equation 



time-scale and the peak luminosity 



ij,46 > 30- 



(63) 



In the specific case of PKS 2155—304, this condi- 
tion implies a lower limit on the jet power, Lj > 
lO'^^MlH.srj'i.sCi erg s-i. Since the Eddington lu- 
minosity has the following value LecI ~ 1.4 x 
10^^AfBH,8 erg s~^, the jet should have a super Edding- 
ton luminosity, unless the bulk Lorentz factor is large, i.e 

In the framework of JRGI scenario, the above condi- 
tions are quite general and do not depend on the radia- 
tion mechanism. Obviously, the available energy has to 
be transferred to gamma-rays by a radiation mechanism 
with an adequate cooling time. It was shown in Section ^ 
that, in general, both leptonic (EIC) and hadronic (pro- 
ton synchrotron) mechanisms can be quite efficient in the 
frameworks of the JRGI scenario. Below we check the 
feasibility of these two mechanisms for the specific case 
of PKS 2155—304; namely, we combine the energy and 
variability constraints with radiation mechanism specific 
limitations (e.g. cooling time). 

4.1. EIC model for PKS 2155-304 

Given the short cooling time in the case of EIC (see 
Equation (^), the required variability can be achieved 
for a relatively low jet luminosity: 



Lj 46 > 0.007 



^"BH,8^ j,1.5 

4/3 2/3 
^2 



(64) 
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On the other hand, the IC scattering should occur in 
the Thomson regime, otherwise the gamma-gamma opac- 
ity would significantly exceed unity (see Equation (pO|)). 
Thus, Equation ( ^8| ) allows the following upper limit on 
the jet luminosity: 



LjA6 < 0.4 



A/2 p6 ,.: 



^11 



(65) 



which appears to be rather close to the lower limit given 
by Equation (|63|). For the sake of clarity, we have com- 
bined all the relevant limitations in Figure 0, where the 
favorable parameters are indicated by the filled region. It 
is expected that the EIC mechanism can be responsible 
for nonthermal emission with the required characteristics 
in the case of sub-Eddington jets. On the other hand, the 
EIC model requires less comfortable bulk Lorentz factor, 
which for the radiation production region should exceed 
Fj > 75. In the framework of the JRGI scenario this 
implies that the star enters into the jet quite far from 
the BH, z > 4 X 10^''AfBH,8 cm. Finally, we note that 
the EIC model does not pose a strict requirement on the 
acceleration process given the relatively small energy of 
the emitting electron. 

4.2. Proton synchrotron model for PKS 2155—304 

One of the most fundamental challenges for proton syn- 
chrotron models in AGN jets is the long cooling time. 
This is often interpreted as an indication of the ext r eme! 



low efficiency of such models. In particular, Bikora (201 



f) 



has compared the synchrotron cooling time to the ex- 
pected adiabatic cooling time in relativistic jets. Based 
on this comparison, it was concluded that this mecha- 
nism was not feasible for nonthermal photon production 
in AGN jets. However, we note that this estimate de- 
pends strongly on the key assumption of the jet magne- 
tization . This param eter has been accepted to be quite 
low by Sikora (2010), whilst in this paper, we rely on 
models of magnetically driven jets, i.e. the magnetiza- 
tion can be very high. Given the large magnetic field, 
proton synchrotron radiation from jets can be quite effi- 
cient in our scenario. 

The cooling time for proton synchrotron radiation. 
Equation (Fsf) , recall for a very luminous jet 



L]A6 > 



500 



8/3 2/3 



8^ j,1.5'/l 



4/3 



(66) 



Note that. Equations (|62|), ( |63| ) and (|66|) together require 
quite an extreme jet luminosity, but the requirement may 
be relaxed in the case of relatively small bulk Lorentz 
factors, i.e., Fj 20. Although even in this case, the 
jet luminosity exceeds the Eddiiigton luminosity. The 
relevantparameter space is shown by the filled region in 
Figure ^ This corresponds to the case when the JRGI 
scenario is realized relatively close to the base of the jet, 
z « 3 X 10^^Mbh,8 cm. Finally, we note that proton 
synchrotron models require a very efficient particle ac- 
celeration with a rate close to the theoretical electrody- 
namic limit, i.e. ry < 10. The corresponding Doppler 
factor can be as high as (5 > 40. One should expect a 
cutoff in the gamma-ray spectrum at around w 1.5?7;j"^ 
TeV. Thus, the observed VHE spectrum, extending up 
to 4 TeV, requires an acceleration parameter of 77 < 10. 



This implies that the proton synchrotron model can be 
realized only in extreme accelerators. 

5. DISCUSSION 
The ultra-short TeV gamma-ray fla res of blazars de 



tecte d in the case of P KS 2155—304 f Aharonian et al 



20071) and Mkr 501 ( [Albert et "al| pOOT^ on 100 s 
timescales represent a serious challenge for current mod- 
els of blazars. This challenge concerns the origin and the 
sites of formation of these flares, the acceleration and 
radiation mechanisms, the hydrodynamics of relativistic 
outflows, amongst others. Since the upper limit on the 
size of production region, of 3 x 10^^r2 cm, is smaller by 
an order of magnitude than the gravitational radius of 
a black hole of mass IO^Mq (which is required to power 
distant blazars), the only way to avoid the situation of 
invoking quite uncomfortable upper limits on the mass of 
the central black hole (as small as 10 ''Mq), is to invoke 
the Doppler boosting. However, this can be realized only 
in the case of an external origin of the processes which 
cause these ultra-short flares. If the flares are initiated 
by disturbances originating from the central black hole 
(e.g., due to internal shocks), the linear size (in the ob- 
server's frame) of the flare production region cannot be 
smaller than the gravitational radius of the black hole, 
independent of the Doppler factor of the jet. The model 
suggested in this work readily solves the problem of con- 
necting the flares to the interactions of the red giant stars 
with the powerful jets. Due to these interactions the 
red giant loses a significant fraction its atmosphere. The 
cloud, accelerated by the magnetically driven jet up to 
a Lorentz factor of F ~ 30, likely separates into many 
small fragments. These "blobs" represent the ideal sites 
for the production of flares, provided that a significant 
fraction of jet energy absorbed by the cloud is converted 
(e.g. due to relativistic shock acceleration or magnetic 
reconnection) to relativistic particles. 

The effective acceleration of particles is a necessary, 
but not a sufficient condition for the interpretation of the 
gamma-ray radiation of blazars. Any successful model 
of TeV blazars require adequate cooling times through 
gamma-radiation; they should be comparable, or often 
even shorter compared to the characteristic timescales of 
other radiative and non-radiative processes. Generally, 
leptonic models of gamma-ray loud blazars, through the 
realization of SSC or external IC scenarios, do provide 
adequate radiation timescales, but at the expense of the 
assumption of a rather weak magnetic field, typically less 
than 1 G, which in powerful blazars (Lj > 10^^ erg s~^), 
is well below the magnetic field in the jet as long as it's 
concerned with sub-parsec distances (see Figure ra). This 
is a quite challenging requirement of the discussion of 
feasibility of, which is generally ignored in the literature. 
In the JRGI scenario suggested here the problem can be 
formally solved assuming that the magnetic field inside 
the blob is much smaller than in the jet. However, in 
the case of the SSC models, this assumption still does 
not allow a relaxation of the second requirement of an 
extremely large jet Lorentz factor, F > 1000. Although 
such Lorentz factors for the bulk motion cannot be ex- 
cludcdQ, in particular at large, > Ipc distances from the 
BH (sec Figure ||) , in the proposed JRGI model it hardly 

^ We should note that it is rather difiicult to reach such a high 



Rapid variability due to Jet-Star interaction 



11 



can work. At such distance the jet ram pressure is not 
suffi cient able to ablate the atmosphere of the star. 



an unavoidable, model-independent conclusion for GRBs 



The requirements on the magnetic field and the jet 
Lorentz factor are more relaxed in the external iC model. 
Nevertheless, one should note that within the JRGI sce- 
nario, the external Compton model has some specific 
features. In order to avoid severe gamma-gamma ab- 
sorption, the Compton scattering should proceed in the 
Thomson regime. This can be fulfilled if the radiation 
region is located at large distances, i.e. regions still with 
quite a large Lorenz factors for the jet, Fj ^ 100. 

One of the main postulates of the JRGI scenario is 
the effective star-jet interaction. This requires the loca- 
tion of the blobs that emit gamma-rays to be at small 
distances from the BH, typically z 
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case of powerful jets, Lj > 
very large magnetic field, B 



lO'^'^ erg s-i 



this implies a 
100 G and a moderate 
Lorenz factor, Fj ^ 20. Both parameters match nicely 
with the interpretation of the TeV gamma-ray flares as 
a result of proton-synchrotron radiation by highly mag- 
netized blobs, formed and accelerated in jet-star inter- 
actions. This model demands the acceleration of pro- 
tons to energies of lO^^eV, and implies the acceleration 
of protons with a rate close to the maximum (theoreti- 
cally possible) rate of, tacc ~ '"l /c This is quite a robust 
requirement, which however, can be provided, in princi- 
ple, by certain acceleration mechanisms. Another chal- 
lenge of the proposed scenario is related to the power of 
the jet. Namely, the proton synchrotron model of TeV 
gamma-rays can be effective, provided that: (i) the mass 
of BH does not significantly exceed M ~ 1O*M0 and; 
(ii) the jet power is not significantly below 10^^ erg s~^. 
An unambiguous implication of these two requirements 
(working in two different directions) is that the jet should 
have a super-Eddington luminosity. Although this could 
seem like quite a dramatic assumption, we note the re- 
quirement of super-Eddington luminosities seems to be 



and also l ikely for powerful gamma-ray blazars (Ghis 



elhni 2011). 



5.1. Stellar density in the vicinity of a SMBH. 

An important question in the suggested scenario is the 
expected rate of the flaring events, which is related to the 
number density of RGs at the relevant jet scales. The jet 
region suitable for the production of the powerful flares 
(similar to the burst detected from PKS 2155—304), can 
be deflned as z < Ipc, and the corresponding side cross 
section of the jet is Sj « ~ lO^^fl-iz^jCm^. Thus, 
the number of flaring event&_per year can be estimated 
as T « S'jV^ib"-- Equation m) provides an estimate for 
the density of RGs required to produce T flaring events 
per year: 



(67) 



Unfortunately, there are no direct measurements of the 
stellar density in the vicinity of BHs. Thus, depend- 
ing on the assumed distribution law, the number of RGs 
in the vicinity of the BH may or may not be sufficient. 
However, we note that studies of possible stellar den- 
sitv profiles in the vicinity of the BH in AGNs (s ee e.g. 
( Bisnovatyi-Kogan et al. 1982 ; Murphy ct al. 1991 ) ) show 
that densities similar to the one required {'^ 10*^ pc~^) 
are rather feasible. Moreover, under the influence of X- 
ray radiation close to BHs, normal stars can drastically 
increase the rate of their stellar wind. Thus, wind-fed 
clouds can be formed during the jet - star interaction. 
This interaction can mimic the interaction of RG atmo- 
sphere with the jet (Dorodnitsyn, private communica- 
tions). Since, the stellar density of normal stars is higher 
up to 2 orders of magnitude than the density of RGs, this 
effect can significantly relax the requirement imposed by 
Equation ( pT\ ) on the stellar density in the vicinity of 
BHs. 



APPENDIX 

DERIVATION OF THE CLOUD DYNAMICS EQUATION IN THE RELATIVISTIC STAGE 

Let us choose the orientation of coordinate systems such that in the observer reference frame, K, the magnetic and 
electric field vectors have the following components B = (0, B, 0) and E = {E, 0, 0), respectively. Since the conductivity 
of plasma is very high, the following condition is held E = VjI3/c, where Vj is jet velocity. The cloud velocity in the 
system K is Vc, and K' is the cloud's momentarily comoving reference frame. Quantities pertaining to the system K' 
are marked by prime. In the system, K' electromagnetic field strengths are 



E' = F, 

B' = rr 



E- —B 
c 



-Fc(uj - Vc)B , 



B-^E]=rr 



1 - 



r.2 



B. 



(Al) 
(A2) 



As it should be, the following ratio holds; E' / B' = where uj is the jet speed relative to the cloud. Energy fiux 
density in system K' and in the laboratory system are related by: 



—E'B' 

477 



(A3) 



Jet ram pressure in the system K' is equal to q' /c^ thus during a differentially small time interval, dt' , the cloud 
momentum increases by a value of dPi, = {q' / c) nr'^ dt' , and the energy increment is a second order value of dt' , i.e 
dEl, ^ 0. In the observer system, one has dE^ = McC^dTc = TcVcdP^, and dt = T^dt' . Thus, the equation of motion 



value of the bulk Lorentz factor, e.g. due to the so called "photon 



breeding mechanism" ( ptern fc Poutanen 2006), which does not 
allow AGN jets with bulk i^orentz tactors exceeding 50. 
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may be expressed as follows: 

Denoting q — Lj/noj^ and assuming Fj ^ 1 and Fc 3> 1, one obtains the equation of motion in the form of Equation ( p^ ) 

(i tA L,rl 

dt \Ti r4 j Auj^cm, ' ^ ' 

THE VARIATION OF 6 IN THE CASE OF BLOB CHAOTIC MOTION 

In the suggested scenario, non-thermal particle acceleration is triggered by RG material blobs caught up in the jet. 
In such a case it is natural to assume that the acceleration sites are closely related to the blobs. Although the particle 
propagation and isotropisation processes may be very complicated, the highest energy particles lose their energy very 
close to the acceleration site. Thus, in what follows we assume that the VHE emission is related to the blobs and the 
radiation boosting factor is determined by the blob velocity. Since the blob velocity can be changed relatively easy, 
sudden changes of the Doppler factor are rather feasible. Hence we have discuss the related changes of the Doppler 
factor and correction function F^.^ i.e. the quantities describing the non-thermal emission intensity, as seen by observer. 

A 

Let us assume that at the moment to, a blob moves relativistically in the direction of the observer with velocity 
V = /3c, thus, the corresponding Lorentz factor is F = l/Vl — ^ 1. Let us consider two reference frames: the 
observer coordinate system K, and the blob comoving system K', where the blob is at rest at to. The z axises are 
selected to be parallel to the systems' relative velocity, V. Let us consider two cases for a change of the blob velocity. 

In the K' reference frame, the blob gains velocity v', such as: 

< = c/3' sin , v'y^O, w;=c/3'cos0. (Bl) 

Here (j> is the angle between the velocity v' and z' axis; /?' = v' /c; and the observer detects the emission radiated along 
z axis. In the laboratory frame, the blob velocity components have the following form: 

c /3'sin0 (3 + l3'coscj) 

Vt — : V,, = V , V, ^ c . B2 

F 1 + cos (/) ' ^ ' 1 + ^/3' cos(/. ^ ' 



Therefore, the blob Lorentz factor is F = 1/ y 1 — (u^ + ij2 -(- vD/c^. Given that the system is in the ultrarelativitic 
regime (in the Appendixes we will assume F ^ 1), one obtains that the Doppler factor is: 



~ 1-1 2Fv/r^ 

F(l - vjc) = Tr> 1 ■ ^^^^ 

i 1 — p' cos 

In the laboratory frame the velocity deflection angle is 50 — (5' sin i/i/F ^ 1. Thus, the Lorentz factor of the radial (i.e. 
along the z axis) motion F^ can be expressed as 



VI - «.7c2 \/l-/?'cOS0 ^ ' 

The radiation intensity is defined by the correction Junction, if the time-scale of the blob velocity variation is longer 
than the typical particle acceleration time. In the case of Fj Fc ^ 1, one obtains 

Tl (1 + /3' cos (/))(! cos 0)3 ' ^ ' 

Otherwise (i.e. if the particle cooling time is long as compared to the time-scale of the velocity change) , the apparent 
luminosity is proportional to the standard factor (5^. Thus, one obtains 



As it can be seen from the comparison of Equations ( pq ) and (Bf), the function r has a weak dependence on the 



velocity change. At the moment to, this function has a value of ro = ^('/')|^/_o ~ 16F . The change of the velocity 
may lead both to an increase and decrease of the function r(0): if = and j3' > 1/3, the correction function value 
is r(0) > 2ro; and if = tt/2 and /?' > 0.54, one has r(7r/2) < ro/2. 
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Another potentially important situation occurs when a blob moves along a helical trajectory, e.g. along the dominant 
magnetic field line. In this case, the additional velocity component is oriented perpendicularly to the averaged velocity, 
i.e. its components may be expressed as the follows: 

< = cP' cos , Vy^ c(3' sin (f> , v', = . (B7) 

Obviously, if an observer is located in the direction of the z' axis, no flux variability is detected. On the other hand, 
if the observer is located slightly off axis, a periodic increase of the flux level is seen. Let us assume that the observer 
is located in the direction of the blob motion for ^ = 0. In this case, the Doppler factor has the following form: 



l-^'2(2cos0-l)- ^^^^ 
The radiation intensity, then, is proportional to the following factor: 

(l-/3'2(2cos0-l))'' 

Obviously, the maximum value of r is achieved when = 0, i.e. when the radiation is emitted towards the observer. 
If the precession velocity fulfills the following condition 



21/4 _ 1 



1/2 



^'> U/^ + l-2cosJ ' 

then at the corresponding moment, the emission intensity can differ by a factor of 2 as compared to the maximum. 
For example, ii cj) = tt/A then the required velocity value is /?' « 0.5. In Figure ^ we show the (/)-angle dependence of 
the normalized intensity for several different values of (3' . 

Another important issue is the time modulation of the emission, as seen in the observer reference frame. Let us 
assume that the azimuthal angle has the following time dependence; (j) = '^i': where w is a constant; and t' is time 
in the system K' . Then a quantity Tq = 27r/a; corresponds to the precession period in the system K' . Obviously, in 
the observer frame the emission intensity should have a different time dependence than the one that can be obtained 
by the substitution of (f> = ujt' into the function shown in Figure ^. Indeed, the radiation emitted at time tr, by the 
source is detected by an observer located ai r = rn at time t, where 

t = tr + \r- r{tr)\/c'^tr + r/c- {nr{tr))/c, (Bll) 



where r{tr) is the emitter position at time tr- Since the constant term r/c can be neglected. Equation (Bll) can be 
reduced to the following form: 

t = tr - {nr{tr))/c. (B12) 

For A'-frame quantities, the azimuthal angle can be represented as = tutr/T. This relationship allows us to express 
t through (j): 

To{l + n 2(3'^ . \ 



(B14) 



By increasing (f> by 27r, the time, t, changes by an amount equivalent to: 

^ _ Toil + 13'^) 
2T 

which is the period of the observed emission. The dependence of the observed intensity, as a function of the observer 
time, is shown in Figure |l^. As it can be seen from this figure, the distribution width decreases with an increase of 
(3'. The emission "half-decay" interval, At/2, is shown in Figure ^as a function of /?'. 
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Fig. 1. — Schematic illustration of the scenario. If a star enters a powerful AGN jet, the outer layers of the star atmosphere are to be 
ablated. Because of the interaction with the jet the ablated cloud expands and gets involved into the jet bulk motion. 
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Fig. 2. — Solutions of Equation ( |l3| ) shown as Fe vs. the relative Lorentz factor of the cloud (left panel) and as Fe vs. the observation time 
(in units to = zo/2DT?c) (right panel). The Lorentz factor of the jet is assumed to be Fj = 30. The foUovfing values of the D-paramoter 
were used: D = 100 (dot-dashed lines), D = 10 (dashed line), D = 1 (solid line) and D = 0.1 (dotted line). 
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z = 10 cm 




Fig. 3. — Sketch of the jet together with characteristic magnetic field strength and bulk Lorentz factors at typical distances from the BH. 
The BH mass and the dimensionless parameter of rotation were assumed to be A/bh = 10**A/q and a Rj 1, respectively. Initially, the jet is 
assumed to be magnetically dominated with the magnetization parameter a > 100. 
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Fig. 4. — Expected SED of blazars in the framework of the proton-synchrotron scenario in the blob's reference frame. The maximum 
at 10^^ Hz is due to synchrotron radiation from secondary (pair-produced) electrons; the maximum at ~ lOOMeV corresponds to the 
synchrotron radiation of primary electrons accelerated in the blob; the maximum at ~ SOOGeV is generated by protons through the 
synchrotron channel. 
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Fig. 5. — Intensive VHE flares (in tlie case EIC) of L.y fa 10^'^ergs ^ and of duration r fa 200 s that are energetically allowed in the 
framework of our JRGI scenario in systems characterized by a jet of luminosity, Lj, powered by a black hole of mass A^bHj are shown by 
the green filled region (we assume Fj = 90 and = 1). The th in s olid line is the limit from Equation (pj]); the thick aajid line is the limit 
from Equation (|65[); the doted line is the limit from Equation (|64[); the dot-dashed line is the limit from Equation ([6^); and the dashed 
line is the Eddingfon limit. 
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Fig. 6. — Intensive VHE flares (in the case EIC) of L.^ ^ lO^'^ergs"^ and of duration r fa 200 s that are energetically allowed in the 
framework of our JRGI scenario in systems characterized by a jet of luminosity, Lj , propagating with bulk Lorentz factor Fj , are shown by 
the green filled region (here we asanme Mbh = 5 x IO^Mq and = 1). The thip_solid line is the limit from Equation (^sj); the thick solid 
line is the limit from Equation (KSh; the doted line is the limit from Equation (f34); the dot-dashed line is the limit from Equation (p2(); 
and the dashed line is the EddingEon luminosity. 
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Fig. 7. — Intensive VHE flares (in the case of proton synchrotron) of L-^ Rs 10'*'' ergs ^ and of duration r f» 200 s that are energetically 
allowed in the framework of the JRGI scenario in systems characterized by a jet of luminosity, Lj, powered by a black hole of mass Mbhj 

; assume Fj = 20 and = 1). The thin solid line is the limit from Equation (k3); the 
; the dot-dashed line is the limit from Equation (p2); and the dashed line is the Edmngton 



which are shown by the pink flUed region 
thick solid line is the limit from Equation (I 
luminosity. 
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Fig. 8. — Intensive VHE flares (in the case of proton synchrotron) of Ri 10*'' ergs~^ and of duration r f» 200 s that are energetically 
allowed in the framework of the JRGI scenario in systems characterized by a jet of luminosity, Lj , propagating with bulk Lorcntz factor Fj , 
are shown by the pink flUed region (we assume AIbh = 5 x IO^Mq and = 1). The thin solidJine is the limit from Equation (p3); the 
thick solid line is the limit from Equation (p9); the dot-dashed line is the limit from Equation (m); and the dashed line is the Edcungton 
luminosity. 
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